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Abstract 

The ratio of branching fractions of the radiative B decays — )■ K*^'y and B^ — ?• cp'y 
has been measured using 0.37 fb~^ of pp cohisions at a centre of mass energy of 
^/s = 7 TeV, cohected by the LHCb experiment. The value obtained is 

-g*°7) _ , „ , ^ ns+0-06 +0.09 
B{BO^^^) -^•^^±"-U«-o.04-o.08> 

where the first uncertainty is statistical, the second systematic and the third is 
associated to the ratio of fragmentation fractions fs/ fd- Using the world average 
for 5(5° K*^-f) = (4.33 ± 0.15) x 10-^ the branching fraction ^(5° (/)-/) is 
measured to be (3.9 zt 0.5) x 10~^, which is the most precise measurement to date. 
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1 Introduction 



In the Standard Model (SM) the decays — )■ K*'^'-f and 5° — > proceed at leading 
order through b s'j one-loop electromagnetic penguin transitions, dominated by a 
virtual intermediate top quark coupling to a boson. Extensions of the SM predict 
additional one-loop contributions that can introduce sizeable effects on the dynamics of 
the transition [T]. 

Radiative decays of the B'^ meson were first observed by the CLEO collaboration 
in 1993 [2] through the decay mode B — )■ K*'y. In 2007 the Belle collaboration re- 
ported the first observation of the analogous decay in the B^ sector, -8° — )■ 07 [3]. 
The current world averages of the branching fractions of B^ — )• K*^'y and B^ — )■ 07 are 
(4.33 ±0.15) X 10-^ and (5.7+?;i) x 10-^ respectively HIS]. These results are m agree- 
ment with the latest SM theoretical predictions from NNLO calculations using SCET 
B{B^^ K*^-f) = (4.3 ± 1.4) X 10-5 and i3(5°-> 07) = (4.3 ± 1.4) x 10-^ which suffer 
from large hadronic uncertainties. The ratio of experimental branching fractions is mea- 
sured to be i3(i?°— )■ i^r*°7)/i3(i?g — )■ 07) = 0.7 ± 0.3, in agreement with the prediction of 
1.0 ±0.2 [6j. 

This paper presents a measurement of B{B^ — )■ K*^'y)/B{B^ — )■ 07) using a strategy 
that ensures the cancellation of most of the systematic uncertainties affecting the mea- 
surement of the individual branching fractions. The measured ratio is used to determine 
B{B^^ (p'j) assuming the world average value of B{B^ ^ K*^^) [1]. 

2 The LHCb detector and dataset 

The LHCb detector [7j is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r] < 5, designed for the study of particles containing 6 or c quarks. The detec- 
tor includes a high precision tracking system consisting of a silicon-strip vertex detector 
(VELO) surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations 
of silicon-strip detectors and straw drift-tubes placed downstream. The combined track- 
ing system has a momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% 
at 100 GeV/c, and an impact parameter (IP) resolution of 20 /xm for tracks with high 
transverse momentum. Charged hadrons are identified using two ring-imaging Cherenkov 
(RICH) detectors. Photon, electron and hadron candidates are identified by a calorime- 
ter system consisting of scintillating-pad and pre-shower detectors, an electromagnetic 
calorimeter and a hadronic calorimeter. Muons are identified by a muon system com- 
posed of alternating layers of iron and multiwire proportional chambers. The trigger 
consists of a hardware stage, based on information from the calorimeter and muon sys- 
tems, followed by a software stage running on a large farm of commercial processors which 
applies a full event reconstruction. 

The data used for this analysis correspond to 0.37 fb"^ of pp collisions collected in the 

^Charge conjugated modes are implicitly included throughout the paper. 
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first half of 2011 at the LHC with a centre of mass energy of ^/s = 7 TeV. ^ K*^'-f and 
— 07 candidates are required to have triggered on the signal photon and vector meson 
daughters, following a definite trigger path. The hardware level must have been triggered 
by an ECAL candidate with > 2.5 GeV. In the software trigger, the events are 
selected when a track is reconstructed with IP > 16, and either > 1.7 GeV/c when 
the photon has > 2.5 GeV or > 1.2 GeV/c when the photon has E't > 4.2 GeV. 
The selected track must form a K*^ or candidate when combined with an additional 
track, and the invariant mass of the combination of the K*^{(j)) candidate and the photon 
candidate is requested to lie within a 1 GeV/c^ window around the nominal B^[B^) mass. 

Large samples of B^ — )■ K*^'-f and B^ — )■ 07 Monte Carlo (MC) simulated events [8] 
are used to optimize the signal selection and to parametrize the B meson invariant mass 
distribution. The pp collisions are generated with Pythia 6.4 [0] and decays of hadronic 
particles are simulated using EvtGen [TU] in which final state radiation is generated 
using Photos [H]. The interaction of the generated particles with the detector and its 
response are simulated using Geant4 [12]. 

3 Event selection 

The selection of both B decays is designed to ensure the cancellation of systematic un- 
certainties in the ratio of their efficiencies. The procedure and requirements are kept as 
similar as possible: the B'^lBg) mesons are reconstructed from a selected ii'*°(0), com- 
posed of oppositely charged kaon-pion (kaon-kaon) pairs, combined with a photon. 

The two tracks from the vector meson daughters are both required to have 
Pt > 500 MeV/ c and to point away from all pp interaction vertices by requiring IP > 25. 
The identification of the kaon and pion tracks is made by applying cuts to the particle 
identification (PID) provided by the RICH system. The PID is based on the comparison 
between two particle hypotheses, and it is represented by the difference in logarithms of the 
likelihoods (DLL) between the two hypotheses. Kaons are required to have DLL^n > 5 
and DLLxp > 2, while pions are required to have DLLi^^r < 0. With these cuts, kaons 
(pions) coming from the studied channels are identified with a ~ 70 (83) % efficiency for 
a ~ 3 (2) % pion (kaon) contamination. 

Two-track combinations are accepted as K*^{(f)) candidates if they form a vertex with 

< 9 and their invariant mass lies within a ±50 (±10) MeV/c^ mass window of the 
nominal K*^{(j)) mass. The resulting vector meson candidate is combined with a photon of 
Et > 2.6 GeV. Neutral and charged electromagnetic clusters in the ECAL are separated 
based on their compatibility with extrapolated tracks [13j while photon and 7r° deposits 
are identified on the basis of the shape of the electromagnetic shower in the ECAL. The 
B candidate invariant mass resolution, dominated by the photon contribution, is about 
100 MeV/c^ for the decays presented in this paper. 

The B candidates are required to have an invariant mass within a ±800 MeV/c^ win- 
dow around the corresponding B hadron mass, to have p^ > 3 GeV/c and to point to a pp 
interaction vertex by requiring IP < 9. The distribution of the helicity angle ^h, defined 
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as the angle between the momentum of either of the daughters of the vector meson (V) 
and the momentum of the B candidate in the rest frame of the vector meson, is expected 
to follow sin^ for B — > V^, and cos^ 9^ for the B — > Vn^ background. Therefore, the he- 
licity structure imposed by the signal decays is exploited to remove B — )■ Vn^ background, 
in which the neutral pion is misidentified as a photon, by requiring that | cos^^hI < 0.8. 
Background coming from partially reconstructed 6-hadron decays is rejected by requiring 
vertex isolation: the of the B vertex must increase by more than half a unit when 
adding any other track in the event. 



4 Determination of the ratio of branching fractions 

The ratio of the branching fractions is calculated from the number of signal candidates in 
the B^^ K*'^j and 4>j channels, 

where N corresponds to the observed number of signal candidates (yield), B((f) — )■ K^K~) 
and B{K*^ — )■ K^n^) arc the visible branching fractions of the vector mesons, fs/fd is 
the ratio of the B^ and -B^ hadronization fractions in pp collisions at y/s = 7 TeV, and 
eBO-^(jrr/^BO-^K*0"f is the ratio of efficiencies for the two decays. This latter ratio is split 
into contributions coming from the acceptance (race), the reconstruction and selection 
requirements (rreco), the PID requirements (rpio), and the trigger requirements (rtrig) : 

= J^acc X Treco X TpiD X Ttrig. (2) 

The PID efficiency ratio is measured from data to be rpiD = 0.787 ± 0.010 (stat), 
by means of a calibration procedure using pure samples of kaons and pions from 
£)*±_^ D°(ir+7r~)7r='= decays selected utilizing purely kinematic criteria. The other ef- 
ficiency ratios have been extracted using simulated events. The acceptance efficiency 
ratio, race = 1-094 ± 0.004 (stat), exceeds unity because of the correlated acceptance of 
the kaons due to the limited phase space in the (f) K^K^ decay. These phase-space 
constraints also cause the (p vertex to have a worse spatial resolution than the K*^ vertex. 
This affects the — > 07 selection efficiency through the IP and vertex isolation cuts 
while the common track cut pt > 500 MeV/c is less efficient on the softer pion from the 
K*^ decay. Both effects almost compensate and the reconstruction and selection efficiency 
ratio is found to be rrcco = 0.949 ±0.006 (stat), where the main systematic uncertainties in 
the numerator and denominator cancel since the kinematic selections are mostly identical 
for both decays. The trigger efficiency ratio rtrig — 1.057±0.008 (stat) has been computed 
taking into account the contributions from the different trigger configurations during the 
data taking period. 

The yields of the two channels arc extracted from a simultaneous unbinned maximum 
likelihood fit to the invariant mass distributions of the data. Signals are described using a 
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Crystal Ball function [1A\ , with the tail parameters fixed to their values extracted from MC 
simulation and the mass difference between the 5° and signals fixed [U]. The width 
of the signal peak is left as a free parameter. Combinatorial background is parametrized 
by an exponential function with a different decay constant for each channel. The results 
of the fit are shown in Fig. [l| The number of events obtained for K*^'y and B^ — )■ 07 

are 1685 ± 52 and 239 ± 19, with a signal over background ratio of 5/5 = 3.1 ± 0.4 and 
3.7 ± 1.3 in a ±3a window, respectively. 




5 I I I I I I I I I I I I I I LJ , , , \ I I , , \ , , , , \ 

4500 5000 5500 6000 5000 5500 6000 

M(KjtY) (MeV/c^) M{KKy) (MeV/c^) 



Figure 1 Result of the fit for the B° K*°-f (left) and B° 07 (right). The black points 
represent the data and the fit result is represented as a solid line. The signal is fitted with a Crystal 
Ball function (light dashed line) and the background is described as an exponential (dark dashed 
line). Below each invariant mass plot, the Poisson residuals [T^ are shown. 



Several potential sources of peaking background have been studied: B^^^ — t- K^tv tt^ 
and B^ — )■ K^K^tt^, where the two photons from the 7r° can be merged into a single 
cluster and misidentified as a single photon, — )■ A*^{Kp) where the proton can be 
misidentified as a pion or a kaon, and the irreducible B^ — )■ K*^^. Their invariant mass 
distributions and selection efficiencies have been evaluated from simulated events and the 
number of predicted background events is determined and subtracted from the signal 
yield. 

B decays in which one of the decay products has not been reconstructed, such as 
S— > {K*^'K^)X , tend to accumulate towards lower values in the invariant mass distribu- 
tion but can contaminate the signal peak. However, their contributions have not been 
included in the fit, and the correction to the fitted signal yield has been quantified by 
means of a statistical study. The mass distribution of the partially reconstructed B de- 
cays is first extracted from a sample of simulated events and the corresponding shape has 
been added to the fit with a free amplitude. The fit is then repeated many times varying 
the shape parameters and the amplitude of the partially reconstructed component within 
their uncertainties. The correction to be applied to the signal yield and its uncertainty at 
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a 95% confidence level are determined from the obtained distribution of the signal yield 
variation. 

The effects of the cross-feed between the two channels, i. e. — > i^*°7 signal misiden- 
tified as — )• 07 and vice-versa, as well as the presence of multiple B candidates per 
event, have also been computed using simulation. The statistical uncertainty due to finite 
MC sample size is taken as the uncertainty in these corrections. 

Table [T] summarizes all the corrections applied to the fitted signal yields, as well as 
the corresponding uncertainties, for each source of background. 

Table 1 Correction factors and corresponding uncertainties affecting the signal yields, in percent, 
induced by peaking backgrounds, partially reconstructed backgrounds, signal cross-feed and multiple 
candidates. The total uncertainty is obtained by summing the individual contributions in quadrature. 



-> ir*°7 5° 07 Ratio 



Contribution Corr. Error Corr. Error Corr. Error 



50 -> K+n-n^ 


-1.3 


±0.4 




< 0.1 


-1.3 


±0.4 




-0.5 


±0.5 




< 0.1 


-0.5 


±0.5 






< 0.1 


-1.3 


±1.3 


±1.3 


±1.3 


A|^^ A*°7 


-0.7 


±0.2 


-0.3 


±0.2 


-0.4 


±0.3 


50 ^ 7^*07 


-0.8 


±0.4 






-0.8 


±0.4 


Partially reconstructed B 


+0.04 


+3.1 
-0.2 


±4.5 


+1.3 
-2.9 


-4.5 


+4.2 
-1.3 


(f)'y/K*^'j cross-feed 


-0.4 


±0.2 




< 0.1 


-0.4 


±0.2 


Multiple candidates 


-0.5 


±0.2 


-0.3 


±0.2 


-0.2 


±0.3 


Total 


-4.2 


+3.2 
-0.9 


±2.6 


+1.9 

-3.2 


-6.8 


+4.5 
-2.0 



The ratio of branching fractions from Eq. [T] is calculated using the fitted yields of the 
signal corrected for the backgrounds, the values of the visible branching fractions [15], 
the LHCb measurement of fs/ fd [13 HE], and the values of the efficiency ratios described 
above. The result is 

— ^ = 1.12 ± 0.08 Stat . 

5 Systematic uncertainties 

The limited size of the MC sample used in the calculation of race, '"reco, and rtrig induces 
a systematic uncertainty in the ratio of branching fractions. In addition, race is affected 
by uncertainties in the hadron reconstruction efficiency, arising from differences in the 
interaction of pions and kaons with the detector and the uncertainties in the description 
of the material of the detector. Differences in the mass window size of the vector mesons. 
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combined with small differences in the position of the K*^{(f)) mass peaks between data 
and MC, produce a systematic uncertainty in rreco which has been evaluated by moving 
the centre of the mass window to the value found in data. The reliability of the simula- 
tion to describe the IPx^ of the tracks and the B vertex isolation has been propagated 
into an uncertainty for rreco- For this, the MC sample has been reweighted to reproduce 
the background-subtracted distributions from data, obtained by applying the sPlot tech- 
nique lin] to separate signal and background components, using the invariant mass of 
the B candidate as the discriminant variable. No further systematic errors are associated 
with the use of MC simulation, since kinematic properties of the decays are known to be 
well modelled. Systematic uncertainties associated with the photon are negligible due to 
the fact that its reconstruction in both decays is identical. 

The systematic uncertainty associated with the PID calibration method has been 
evaluated using MC simulation. The statistical error due to the size of the kaon and pion 
calibration samples has also been propagated to rpiD. 

The systematic effect introduced by applying a B mass window cut of ±800 MeV/c^ 
has been evaluated by repeating the fit procedure with a tighter B mass window reduced 
to ±600 MeV/c^. 

Table |2] summarizes all sources of systematic uncertainty, including the background 
contributions detailed in Table [T} The uncertainty on the ratio of efficiency-corrected 
yields is obtained by combining the individual sources in quadrature. The uncertainty on 
the ratio fs/ fd is given as a separate source of uncertainty. 

Table 2 Summary of contributions to the relative systematic uncertainty on the ratio of branching 
fractions. Note that fs/ fd is quoted as a separate systematic uncertainty. 



Source 


Uncertainty 


(%) 


Acceptance (race) 


±0.3 




Selection (rreco) 


±1.4 




PID efficiencies (rpm) 


±2.7 




Trigger (rtrig) 


±0.8 




B mass window 


±0.9 




Background 


+4.5 
-2.0 




Visible fraction of vector mesons 


±1.0 




Quadratic sum of above 


+5.4 
-3.3 




fs/ fd 


+7.9 
-7.5 





Besides fs/ fd-, the dominant source of systematic uncertainty is the imperfect mod- 
elling of the backgrounds due to partially reconstructed B decays. This specific uncer- 
tainty is expected to be reduced when more data are available. 



6 



6 Results and conclusions 

In 0.37 fb^^ of pp collisions at a centre of mass energy of ^/s = 7 TeV the ratio of branching 
fractions of — )■ K*^'j and — ?■ 07 decays has been measured to be 

Tl^^-^^^y ^ ^-^^ ^ ^o:S8(/.//.) 

in good agreement with the theoretical prediction of 1.0 ± 0.2 [6]. 
Using ^ ir*°7) = (4.33 ± 0.15) x 10"^ [1], one obtains 

07) = (3.9 ± 0.5) X 10~^ 

(statistical and systematic errors combined), which agrees with the previous experimental 
value. This is the most precise measurement of the B^ — )■ 07 branching fraction to date. 



Acknowledgements 

We express our gratitude to our colleagues in the CERN accelerator departments for 
the excellent performance of the LHC. We thank the technical and administrative staff at 
CERN and at the LHCb institutes, and acknowledge support from the National Agencies: 
CAPES, CNPq, FAFERJ and FINER (Brazil); CERN; NSFC (China); CNRS/IN2F3 
(France); BMBF, DFG, HGF and MFC (Germany); SFI (Ireland); INFN (Italy); FOM 
and NWO (The Netherlands); SCSR (Foland); ANCS (Romania); MinES of Russia and 
Rosatom (Russia); MICINN, XuntaGal and GENCAT (Spain); SNSF and SER (Switzer- 
land); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We also acknowl- 
edge the support received from the ERG under FF7 and the Region Auvergne. 



References 

[1] S. Descotes-Genon, D. Ghosh, J. Matias, and M. Ramon, Exploring new physics in 
the C7-C7' plane, JHEF 6 (2011 ) 099, arXiv: 1104. 3342; T. Gershon and A. Soni, 
Null tests of the Standard Model at an international Super B factory, |J. Phys. G33 
[(2007) 479, arXiv:hep-ph/0607230 ; F. Mahmoudi and M. R. Ahmady, Constraints on 
mSUGRA from isospin asymmetry in B ^ K*'y, AIP Conf. Proc. 903 (2007) 283, 
^axXiv:hep-ph/0610144; W. Altmannshofer, P. Paradisi, and D. M. Straub, Model- 
independent constraints on new physics in b ^ s transitions, arXiv:1111.1257 

[2] CLEO collaboration, R. Ammar et ai. Evidence for penguin- diagram decays: first 
observation of B ^ iT* (892)7, Phys. Rev. Lett. 71 (1993) 674, 

[3] Belle collaboration, J. Wicht et al.. Observation of B^ — >■ 07 and search for B^ — )• 77 
decays at Belle, Phys. Rev. Lett. 100 (200 8) 121801, ,arXiv:0712.2659, 



7 



[4] 
[5] 



Heavy Flavor Averaging Group, D. Asner et ai, Averages of h-hadron, c-hadron, and 
T-lepton Properties, arXiv:1010.1589 . 

BaBar collaboration, B. Aubert et ai, Measurement of branching fractions and CP 
and isospin asymmetries in B K*{892)'y decays, Phys. Rev . Lett. 103 (2009)] 
211802, axXiv:0906. 21771 Belle collaboration, M. Nakao et ai, Measurement of the 



B K*^ branching fractions and asymmetries, Phys. Rev. D69 (2004 ) 112001| 
CLEO collaboration, T. E. Coan et ai. Study of exclusive radiative B meson decays, 
|Phys. Rev. Lett. 84 (2000) 5283, 



[6] 

[7] 



A. Ali, B. D. Pecjak, and C. Greub, B ^ V-f decays at NNLO in SCET, Eur. Phys. 



J. C55 (2008) 577, arXiv:0709.4422 



LHCb collaboration, A. Augusto Alves Jr et ai. The LHCb detector at the LHC, 



JINS'T 3 (2008) S08005 



M. Clemencic et ai. The LHCb simulation application, Gauss: design, evolution and 
experience, [Journal of Physics: Conference Series 331 (2011) 032023| 



[9] 
[10] 

[11] 
[12] 

[13] 

[14] 

[15] 
[16] 
[17] 



T. Sjostrand, S. Mrenna, and P. Skands, PYTHIA 6.4 physics and manual, JHEP 
|05 (2006) 026, arXiv:hep-ph/0603175 . 



D. J. Lange, The EvtCen particle decay simulation package, Nucl. Instrum. Meth. 
pV462 (2001)152] 



E. Barberio and Z. Was, PHOTOS - A universal Monte Carlo for QED radiative 
corrections: version 2.0, Comp ut. Phys. Commun. 79 (1994) 291^ 



S. Agostinelli et ai, CEANT4 — a simulation toolkit, Nucl. Instrum. Meth. A506 



(2003) 250^ 



O. Deschamps, F. Machefert, M.-H. Schune, G. Pakhlova, and I. Belyaev, Photon 



and neutral pion reconstruction, LHCb-2003-091 



T. Skwarnicki, A study of the radiative cascade transitions between the Upsilon-prime 
and Upsilon resonances. PhD thesis, Cracow Institute of Nuclear Physics, 1986. 
DESY-F31-86-02. 



Particle Data Group, K. Nakamura et ai. Review of particle physics, J. Phys. G37| 
|(20l0) 075021| 

S. Baker and R. D. Cousins, Clarification of the use of chi-square and likelihood 



functions in fits to histograms, Nucl. Instrum. Meth. A221 (1984) 437. 



LHCb Collaboration, R. Aaij et ai. Determination of fs/fd for 7 TeV pp collisions 
and a measurement of the branching fraction of the decay Bd — )■ D^K^, Phys. Rev.] 
Lett. 107 (2011) 21180H|arXiv:1106.4435| 



8 



[18] LHCb collaboration, R. Aaij et ai, Measurement of b hadron production fractions in 



7 TeVpp collisions, Phys. Rev. D85 (2012) 032008, arXiv: 11 11.2357 



[19] M. Pivk and F. R. Le Diberder, sPlot: A statistical tool to unfold data distributions, 
,Nucl. Instrum. Meth. A555 (2005) 356, .arXiv:pliysics/0402083, 



9 



